Polymer dynamics in chaotic flows with strong shear component 



in 
o 
o 

(N 
G 



U 
d 



> 
in 

o 

o 
in 
o 



> 

x 



K. S. Turitsyn 

Landau institute for Theoretical Physics, Moscow, Kosygina 2, 119334, Russia 
Theoretical Division, LANL, Los Alamos, NM 87545, USA 
(Dated: February 8, 2008) 

We consider the internal dynamics of the polymer molecule which is injected in the chaotic flow 
with strong mean shear component. The flow geometry corresponds to the recent experiments on 
the elastic turbulence (Groisman, Steinberg 2000). The passive polymer in such flows experiences 
aperiodic tumbling. We present a detailed study of the statistical properties of such polymer dy- 
namics. First we obtain the stationary probability distribution function of the polymer orientation. 
Secondly we find the distribution of the time periods between consequent events of tumbling, and 
finally we find the tails of the polymer size distribution. 
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Hydrodynamics and rheology of dilute polymer solu- 
tions has attracted much theoretical and experimental 
attention recently. Addition of small amount of poly- 
mers to ordinary liquid leads to crucial changes of liquid 
properties. One of the most famous effects of this type 
is the phenomenon of drag reduction. The addition of 
few parts per million of long-chain polymer molecules 
produces a dramatic reduction of the friction drag. Al- 
though this effect was first observed by Toms in 1949 
Q, there is still no rigorous theory, explaining the phe- 
nomena. The qualitative description was proposed by 
Lumley 0, @], but no quantitative theory is available. 
Another spectacular phenomena, observed in dilute poly- 
mer solutions is the effect of elastic turbulence, discov- 
ered recently by Groisman and Steinberg 0, Ej ■ In this 
experiment a chaotic fluid motion was observed in the 
system with small Reynolds number Re <C 1. Obviously, 
such behavior can not be observed in Newtonian liquids, 
where the flow should be laminar, so the chaotic flow 
is generated by the elastic instabilities of polymer so- 
lution. The dynamics of polymers and possible mecha- 
nisms, explaining the chaotic state were studied in the 
recent theoretical works 0, 0, @ • It was proposed that 
elastic instabilities occur because of elongation of a sin- 
gle polymer in external flows. The analysis of the single 
polymer dynamics in chaotic flows, shows that the transi- 
tion between two qualitatively different types of behavior 
of polymers can be observed in such system. This transi- 
tion is called coil-stretch transition. It separates the dy- 
namics in weak flows, where polymer molecules remain in 
coiled state most of the time, and strong flows, where the 
molecules become substantially elongated. The measure 
of the flow strength is given by the Weissenberg number 
which is the product of the characteristic velocity gra- 
dient and the polymer relaxation time. More precisely 
the Weissenberg number is defined as Wi = At, where 
A is the largest Lyapunov exponent, associated with the 
flow, and r is the relaxation time of the slowest poly- 
mer excitation mode. The coil-stretch transition occurs 
at Wi = 1. With the development of novel optical meth- 
ods a number of high quality experimental observations 
focusing on resolving dynamics of individual polymers 



(DNA molecules) subjected to a non-homogeneous flow 
have been reported [3II1II1II3- 

It made possible the 
direct observation of coil-stretch transition [13j . However 
in all theoretical works, describing the polymer dynam- 
ics in chaotic flows one of the important assumptions was 
the isotropic statistics of the velocity field. At the same 
time in the experiments on clastic turbulence the flow 
consists of regular (shear- like) and chaotic components, 
the latter is relatively weak in comparison with the for- 
mer one. Thus, it became an important task to generalize 
the known theoretical results on the flows of such type. 

In this paper we study the dynamics of the DNA poly- 
mers injected in the chaotic flow. The flow consists of 
large stationary shear component with the small chaotic 
part. It is assumed that the concentration of the in- 
jected polymers is small enough so that their back re- 
action on the flow is negligible. It was shown in 0] 
that the polymer spends most of the time stretched in 
the shear direction and experiences aperiodic tumbling 
towards the opposite direction. The main aim of this 
paper is to present a detailed investigation of the statis- 
tical properties of such polymer dynamics and to make 
some predictions which can be checked experimentally. 
First, we obtain the exact expression of the polymer 
direction distribution. We show that the body of the 
angular distribution function is located in the region of 
small angles, which correspond to the polymer stretched 
in the shear direction. However, the tails of the PDF 
are algebraic, so the fluctuations of the polymer direc- 
tion are very strong. Secondly we study the statistics of 
the time periods between consequent events of polymer 
tumbling. It is shown that this PDF has an exponential 
tail at large tumbling times, and two different asymp- 
totic regimes in the region of very small times. Finally, 
for polymers, which are below the coil-stretch transition 
we obtain the asymptotic of the polymer size distribu- 
tion function, which is also algebraic. We also show, that 
mean-shear component leads to a significant broadening 
of the polymer size distribution. This effect is surpris- 
ing at the first sight, because the shear component itself 
doesn't lead to an exponential polymer growth, and thus 
can not lead to the algebraic tails of the polymer size 
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distribution. However the combined effect of the strong 
regular shear and chaotic velocity components leads to a 
significant enhancement of the polymer growth rate. 

The overall plan of this article is following: in first 
section we describe the experimental setup in which the 
elastic turbulence is observed, then we discuss the as- 
sumptions of the models, describing internal polymer dy- 
namics and flow statistics. In the main sections we give 
the detailed analysis of the polymer dynamics, including 
the tumbling phenomenon j 1 31 ] and statistical properties 
of polymer elongation. While in this paper we present 
only analytical predictions, some of them have been suc- 
cessfully checked numerically in other paper 



I. EXPERIMENTAL SETUP 

The classical effect of well-developed turbulence is ob- 
served in the systems with large values of Reynolds num- 
ber Re = VL/v ^S> 1, where V is the characteristic value 
of the flow velocity, L is the characteristic system size, 
and v is kinematic viscosity. In this case the chaotic 
motion is generated by the large nonlinear terms in the 
Navier-Stokes equation ^(|. In the recently discovered 
elastic turbulence phenomena the chaotic behaviour is 
observed at extremely small Reynolds numbers Re <C 1, 
and is due to elastic instabilities of the polymer solution. 
The experimental setup in which this effect has been ob- 
served looks as following: a swirling flow is generated 
between two coaxial disks of the radius R and with the 
gap d <C R. The upper disk is being rotated with the an- 
gular frequency f2, while the lower one is stationary. For 
small frequencies or usual Newtonian liquid the following 
laminar stationary flow will be generated: 

V4, = nrz/d, V r = V z = 0, (1) 

Here 4>,z,r are the usual cylindrical coordinates with the 
center coinciding with the center of the lower disk. In 
the dilute polymer solution a small chaotic component of 
the velocity field v appears after some critical value of fl 
(but still at small Reynolds numbers). Currently existing 
theoretical works are unable to describe the statistical 
properties of such flow. Thus, it is of interest to study 
this properties experimentally. One of the experimental 
possibilities of such study is the observation of dynamics 
of single long polymer molecule (such as DNA) . It should 
be stressed, that the observed molecules are not related 
to the polymers, which are dissolved in the liquid, and 
which generate the chaotic flow. The concentration of the 
observed DNA molecules is negligibly small, so they can 
be treated as passive objects. In the following text we 
won't analyze the origins of the elastic turbulence, and 
will denote by the polymers the passive observed DNA 
molecules. 
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FIG. 1: Schematic figure explaining polymer orientation ge- 
ometry. 

II. POLYMER AND CHAOTIC FLOW MODELS 

For typical polymers the inertial effects can be ne- 
glected, and one can assume that the polymer simply 
follows the lagrangian trajectories of the velocity field. 
In this case the velocity field affects the internal polymer 
dynamics in two ways: firstly the local velocity gradients 
lead to the polymer stretching and rotation, secondly, if 
the flow parameters are spatially inhomogeneous (which 
happens for example in the experimental setup described 
above), the lagrangian advection can lead to the time 
dependent local velocity gradients, corresponding to the 
stationary velocity component. In our analysis we ne- 
glect the second effect. In the case of the elastic turbu- 
lence experimental setup described above, the advection 
in the radial direction, which leads to the variations of 
the local shear rate is determined only by the chaotic 
component which is assumed to be small. It means the 
that significant changes of the local shear rate, due to la- 
grangian advection occur on the time scales, much larger 
than the characteristic time scales associated with the 
internal polymer dynamics. Therefore in the adiabatic 
approximation one can assume the local shear rates, cor- 
reposding to the regular velocity component to be con- 
stant. Polymer extension can be described by the simple 
dumb-bell model, where the end-to-end separation vector 
R satisfies the following equation 0, : 

d t Ri = RjVjVi - j(R)Ri + Q (2) 

where the relaxation rate 7 is the function of the abso- 
lute value R of the vector R, the velocity gradient VjVi 
is taken at the molecule position, and Q is the thermal 
Langevin force. We assume that velocity field is large- 
scale, so it can be assumed to be smooth on the polymer 
size. This assumption justifies, the linear approximation 
of the velocity field used in J5J. In this paper we discuss 
two different situations. Firstly when the Lyapunov ex- 
ponent, corresponding to the velocity field is larger than 
the polymer relaxation rate, so the coil-stretch transition 
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has already occurred, the polymer becomes nonlinearly 
stretched so, one can neglect the thermal forces Q act- 
ing on him 0| . One can introduce the polymer direction 
vector rii — Ri/R, which is in this case described by the 
following equation: 

dtrii = n j {8 ik - nin k )ViV k (3) 

One can sec that the direction evolution is completely 
decoupled from the dynamics of the polymer size R. In 
the second situation beyond the coil-stretch transition, 
one can not neglect the thermal force d and the dynamics 
of the polymer direction becomes more complicated. We 
will study only the statistics of the polymer size R in this 
case. If the average polymer size is much smaller than 
it's nonlinear length, one can assume the relaxation "f(R) 
to be constant, in which case the equation becomes 
constant and can be easily analyzed. 

Next, it is important to discuss how the chaotic veloc- 
ity component is modeled. The statistical properties of 
the velocity field, observed in the chaotic turbulence are 
not known neither from the experimental nor theoretical 
points of view. In this paper we will study the simplest 
model when the velocity field consists of the strong sta- 
tionary shear component and of the weak short correlated 
chaotic component. Under this assumptions the velocity 
gradient matrix has the following form: 

VjVi = sS iX SjY + (4) 
((Ty^K^O) = DS(t - fO(4Mii - S U 8 kj - SySu) (5) 

Here we assume that the shear flow occurs in the XY 
plane (note that this plane is stationary in the rotating 
frame, associated with the polymer). The exact form 
of the correlation function assumes the isotropy of 
velocity component, however as we will see this assump- 
tion is not important, because in the case of strong shear 
component s D the polymer spends most of the time 
stretched in the X directions, and it's angular dynamics 
is determined only by the Y component of the chaotic 
velocity field. 

In order to simplify the equations on the polymer di- 
rection we parameterize the vector n as shown in the Fig. 
^ The Eq. © acquires the following form 

d t <t> = -ssin 2 4> + £ , (6) 
^ = - S ^sin0 + e e , (7) 

where £^ and t;g are zero mean random variables related 
to chaotic components of the velocity gradient. The 
statistics of both £^ and can be obtained from the 
correlation function JSJ: 

m)( e (t'))=4,DS(t-t') (8) 

mW)) = ^rAt-t') (9) 



III. STATISTICS OF POLYMER DIRECTION 
A. (^-angle distribution 

It follows from © that for the stretched polymers 
the angular dynamics is decoupled from the dynamics 
of the polymer size. In the case of s ^> D the regular 
component will play the main role in the evolution of 
the polymer direction. For vanishing chaotic component 
D = the deterministic polymer dynamics can be eas- 
ily analyzed: there are two semistable equilibrium states 
</>i,2 = 0, 7r; #1,2 = 0, and the polymer direction vector 
n asymptotically approaches one of these points depend- 
ing on it's initial orientation. However when the angle 
between polymer and equilibrium direction becomes suf- 
ficiently small, the chaotic components t; k can not be 
neglected, and the polymer dynamics becomes stochas- 
tic. After some time the chaotic component transfers the 
polymer into instable region, and regular velocity quickly 
(on times of order s _1 ) transfers it to the opposite equi- 
librium direction. Due to stochastic nature of the chaotic 
velocity component one will observe random a-periodic 
tumbling of polymer. This phenomena was qualitatively 
analyzed in [14| for general velocity statistics. In this 
work we focus on the situation where the chaotic flow is 
short correlated, such that it's characteristic correlation 
time t v is small compared to the time scale, associated 
with the tumbling effect which as we will show can be es- 
timated as T t = (Z?s 2 ) -1 / 3 ^> s _1 . This assumption will 
allow us to obtain some rigorous results on the stationary 
and dynamical statistics of polymer. 

When the chaotic component is weak enough D <C s 
the polymer spends most of the time in the stochastic 
regime, close to the equilibrium point, so that it's angles 
are small 9, <f> <C 1 (we will analyze only one equilibrium 
point, because of the symmetry n — * — n). In this case 
one can set 9 = in the correlation function @, so the 
dynamics of angle <p becomes decoupled from everything 
else, and one can write the corresponding Fokker-Planck 
equation: 

[d t - sd^ sin 2 (f> - 2DdJ\ P(t, 0) = (10) 

One of the important questions that should be discussed 
here are the boundary conditions, that should be used 
with this equation. The equation is invariant under the 
transformations (j> — ► <j> ± tt. Therefore it would be nat- 
ural to use periodic boundary conditions P(t,—w/2) = 
P(t,%/2). In this case there will exist an asymptotic 
stationary solution of (|10|l P st (</>). Obviously, all angles 
differing to the value divisible by 7r are identical to each 
other in this solution. Another possibility is to use tradi- 
tional boundary conditions V(t, 00) = V(t, —00) = 0. In 
this case the angles 4> and 4> + irk are not equivalent, and 
one can study the statistics of the number of polymer 
rotations. The main disadvantage of such approach is 
that there is no stationary solution of the Fokker-Planck 
equation, because the PDF is widening and drifting con- 
stantly. However both of the approaches lead to the same 
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physical results, the following identity holds between two 
different PDFs: 



P(t,0=£>(t,^ + 7rA;) 



(11) 



In order to find the stationary PDF P s t{6) we rewrite 
the Fokker-Planck equation in the following form: 



d <t ,U- 1 (<f>)d <t ,U{<f>)P st {<f>) = Q 
U{4>) — exp 



s s 

6 sin 26 

4ZT 8D V . 



(12) 
(13) 



Trivial integration leads to the following expression for 
PDF: 



Pst(<P) 



D 



>exp 



s 

AD 



sin</>cos(0 — 2<f))) 



(14) 

where ui is the average rotation frequency of the poly- 
mer, which is determined from the normalization condi- 
tion J Q P st {<f))d(j) — 1 and is given by 



£>exp £ 



8D, 



TT 2 I ix (x)I- ix (x) ' 



(15) 



where x — s/8D. In the most interesting case s ^> D, the 
PDF will be localized at small angles 6 ~ (D/s) 1 / 3 <C 1 
and all expressions are significantly simplified: 



2\l/3 



4 • 3 1 /6r(7/6) v ^F 



Pst(4>) 



D 



>exp 



— 6(6 - 26) 2 - 

8D 1 ' 24L> 



(16) 



(17) 



One can see, that the PDF is asymmetric and has wide 
algebraic tails: 



P s t{6) 



1 



1/3 



^3 



1/3 



16 • 3 1 /6r(7/6)V^ : V-D 



r(i/6) 

/ S \ 2 /3 1 



(18) 
(19) 



where the last asymptotic is valid in the intermediate 
region {D/s) 1 / 3 -C \6\ <C 1 The asymptotic P st cx 
corresponds to the non-zero probability flux at \6\ — > oo 
which is a manifestation of the described above tumbling 
effect. Positive value of the average angle 6 shows, that 
the polymer spends most of the time in the region > 
in agreement with the qualitative analysis presented in 
the first sections of this paper. One can see, that the 
polymer spends most of the time in the region of small 
angles, therefore the only relevant velocity component v z . 
The assumption of isotropic statistics of chaotic velocity 
component is not therefore significant for the qualitative 
results of this paper. 



B. Tumbling time statistics 

In this section we will calculate the probability distri- 
bution function of the time intervals between consequent 
tumblings. Such PDF can be directly measured exper- 
imentally. For this problem it is naturally to use the 
non-stationary PDF V(t, 6). We will define the tumbling 
process by a polymer direction "trajectory" starting at 
6 = tt/2 and reaching 6 — —ir/2 at time T. In this case 
the probability of finding the polymer inside this region 
is given by 



P(t) 



tt/2 



tt/2 



(20) 



Where the initial condition is V(t, 6) — 5(6 — tt/2 + 0), 
sop(0) = 1. We substitute V(t,6) = U 1/2 {6)^{t,6)- In 
this case the evolution of iff is determined by the one- 
dimensional Schroedinger equation in imaginary time: 



ft* 



H = —2Dd\ H sin 4 6 — s sin 6 cos < 

* 8L> 



(21) 
(22) 



Now it is possible to use the quantum-mechanical anal- 
ogy. The hamiltonian H formally describes a particle in 
periodic potential with the period n. General solution of 
this problem looks like 

*(t,6)=J2 I d P ^ np (6)K P ^/2)exp(-E n (p)t), (23) 

n 

where p is particle's quasimomcntum, and index n enu- 
merates the Brillouin zone number. In this potential the 
classical minimums are separated by large barriers. And 
for s D one can use the tight-binding method [T^ |. 
The following approximate relations hold for the spec- 
trum: 



En(p) 



v cos(-7rp) 



®np(<t>) = ^2exp(i7rkp)ip n (6 



kit) 



(24) 
(25) 



where ip n , e n is the spectrum, which is formed near clas- 
sical minimums, when the tunneling processes are ne- 
glected, v is exponentially small band width. Therefore 
at large times the main asymptotic of p(t) will be deter- 
mined by the ground state energy £q : 



p(t) oc exp(— eot), t —>■ oo 



(26) 



One can easily check that this energy is given by eo = 
c(Ds 2 ) 1/ ' 3 , where c is a constant of order unity. Indeed, 
the classical minima is situated in the region of small an- 
gles 1 0| <C 1, so one can use Taylor expansions of trigono- 
metric functions. After the substitution 6 = (D/s) 1 / 3 
one obtains the Hamiltonian 



H = (Ds z ) 



2x1/3 



(27) 
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The operator in square brackets contains no dimension- 
less parameters, and therefore it's eigenvalues will be of 
order unity. The body of the PDF is also situated in the 
region of tumbling periods of order T ~ (Ds 2 ) -1 / 3 . Left 
tail of tumbling time PDF T <C (Ds 2 ) -1 / 3 is determined 
by rare trajectories, which turn the polymer by angle 7r at 
small times T . In order to find the optimal form of such 
trajectories we will functional integral representation of 
the transition probability 



p(T) oc J T>(f>exp 



1 

~8D 



(28) 



The integration is taken over trajectories with bound- 
ary conditions 0(0) = 7r/2, (f>(T) = —tt/2. For small 



T <C (Ds 2 ) -1 / 3 the probability is determined by the ac- 
tion A on the optimal trajectory p(T) ~ exp(— A) with 
exponential accuracy. Variation of the effective action 
leads to the following equation on instanton: 



s 2 sin 3 i 



(29) 



This mechanical problem can be easily solved, and one 
obtains the following relation between the tumbling time 
and effective particle energy: 



T 



tt/2 



-tt/2 \J2E + s 2 sin 4 . 



E(2E + s 2 ) 



1/4 



K 



E 



AE + 2s 2 



where K(x) is elliptical integral of the first kind. The action in this case is given by 



A 



ET 
1~D 



s 

AD 



d(f> 



V2E~- 



ET 

w 



3ns z 



32DV2E 



" 1 1 2' 4' 4' 2' ' 21: 



(30) 



(31) 



r 



Here we omitted the constant term 2 s J <ft sin 2 <pdt = 
7ts/(8D) because it is cancelled by the normalization con- 
stant. Because of additional time scale s _1 there are two 
different asymptotes of pit) . In the case of sT <C 1 one 
has s < \p2E and E = n 2 /(2T 2 ). In this case 



A = 



7T 



(32) 



8DT 

1 < T < (Ds 2 )- 1 / 3 the en- 



In another limiting case s 

ergy is given by E = 8K 4 (l/2)/(s 2 T 4 ) and the action 
has the following form: 



A = 



2K 4 (l/2) 



(33) 



3L>s 2 T 3 

The intermediate asymptotic l|33|l is a function of the 
product Ds 2 T 3 because it is determined by the dynam- 
ics in the region of small angles. Therefore it will be 
universal, in the sense mentioned above. In contrary, the 
asymptotic at the smallest times l|32ll does not depend on 
s at all, because such small times can be reached only due 
to very rare fluctuation of chaotic velocity field. There- 
fore, this asymptotic strongly depends on the assumption 
of isotropic velocity statistics, and is not universal. 



(D/s) 1 / 3 -C <C 1. The first one is coming from the re- 
gion of deterministic regions where <j> ~ 1, where the an- 
gular dynamics is determined by the regular terms. The 
algebraic tail in this case is proportional to 6~ 2 , however 
there is also a non-universal algebraic part, which comes 
from the stochastic region of ~ (D/s) 1 ^ 3 and is deter- 
mined by the statistical properties of the random velocity 
field. In this section we will be interested in this part, 
and will obtain the connection between the exponent and 
the entropy function of the random velocity process. In 
the region 6 <C 1 the equation Q can be easily solved: 



9(t) 



dTexp 



sin20(£')dt' &(*-t) 



. (34) 

As we already know the random process <f>{t) is stationary 
and independent from £g(i). In this case the expression 
for 6 can be rewritten in the form 



C. #-angle distribution 

As it was shown in there are two contributions 
to the intermediate right tail of the 0-angle distribution 



drexp(-g(r))^e(r) 



B^) = \f sin20(i)di 



(35) 
(36) 
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In order to obtain the PDF P{0) we first average over 
the noise 



2A 



P(6\q) = (27rA)~ 1/2 exp 

A = AD / dTexp[-2g(r 
Jo 



(37) 
(38) 



Here P{0\g) is the PDF of 9 for a fixed realization of the 
process g(t). Due to positive value of (g) ~ (Ds 2 ) 1 / 3 
the body of P(0) is situated in the region of small an- 
gles 9 ~ (D/s) 1 / 3 < 1. Tails of the PDF are deter- 
mined by large deviations of negative g(t). Assuming 
that the process g(t) reaches it's most negative value to 
the time t* , such that g(r*) = —g* and p* ^> 1 one 
can estimate the value of A with exponential accuracy 
as A ~ (D/s) 1 / 3 exp(2g*). The characteristic correlation 
time of g(t) is r c = (Z3s 2 ) -1 / 3 , so for large r* 3> r c one 
can use the results of large deviations theory [2£j] which 
predict the following scaling for the tails of p* PDF: 



P(p*|T*) = exp 



T n T* 



(39) 



where S(x) is the entropy function, which form can not 
be found analytically in general case. One can now find 
the most probable time r* by maximizing the above 
probability over r*. This leads one to the expression 
t* = r c g* jx* , where x is found from the following equa- 
tion: 



S(x*) = x*S'{x*) 



(40) 



The entropy function is of order unity, so one can expect 
the same from x*. The asymptotic of g* PDF is therefore 
given by 



P(g*) ~ exp(- e *SV)) 



(41) 



and after averaging (|37|l over g* one obtains the following 
asymptotic of 9 PDF: 

P{9) cx \6\- sl{ - x '\ (D/s) 1 ' 3 « \0\ « 1 (42) 

One can see that the tails are algebraic as in the case of 
the angle </>, however the exponent is now non-universal, 
and depends on the statistical properties of the velocity 
field. 



IV. STATISTICS OF POLYMER ELONGATION 

The tails of the polymer size PDF can be studied in the 
similar way as in section All CI The explicit solution of 
the dynamical equation [3 in the case of linear relaxation 
force has the following form 

/>OC 

Ri(t)= dtexp[-7(t-t')]^(i,00(i') (43) 



W = Texp 



(44) 



where cry = Vj«, is the velocity gradient matrix. In 
order to obtain polymer size PDF we first average over 
the thermal Langevin force £i(t); 



P(R\a) oc exp 



^R^R 

2 



(45) 



I — K I dt'W 1 (t')W(t) exp(-2 7 i) (46) 
Jo 

where W(t) = W(t,0), and P(R\a) stands for the PDF 
with a fixed realization of the process a(t). At large 
enough times t ^> r c the eigenvalues of the matrix W T W 
become widely separated and the absolute value of the 
end-to-end vector R is determined by the largest eigen- 
value I\\ 



P{R\d) oc exp 



2h 



(47) 



It can be easily shown (see e.g. |21|) that for large times, 
when the eigenvalues of W T W are widely separated 
Ai ^> A2 3> A3 the dynamics of the largest one Ai = 
exp(2p) is described by the following equation: 



p = - cos 2 9 sin 2(j) + 6D + £ p , 



(48) 

(49) 
(50) 



where £ p is obtained from the chaotic velocity correla- 
tion function ©. The eigenvalue l\ is then given by the 
expression 



h = k / dt exp [2p{t] - 2jt] 



(51) 



Like in the previous section p(t) is an integral of station- 
ary random process with the correlation time of order 
t c = (ZXs 2 ) -1 / 3 , and large deviations of I\ are deter- 
mined by the large deviations of pit). Assuming that 
the integral H51J) is determined by one saddle point r* 
and can be estimated like I\ cx exp(2p* — 2jt*) where 
p* = p(r*). The asymptotic of p* PDF for a fixed value 
of t* has the following form: 



P(p*|T*) oc exp 



T -s n 



P T c 



(52) 



One can now find the optimal value of r* = r c p* jx. The 
coefficient x is found from the following equation: 



Spix) - xS'Jx) + JTcS'Jx) 



(53) 



As long as we study the linear region beneath the coil- 
stretch transition we have 7T C > 1. The tail of PDF 
will be algebraic like in case of 9 angle: P{R) oc R^ 1 ^", 
and the value of a can be determined for large values of 
7T C 3> 1. Large deviations of R are determined by the 
region where thermal Langevin forces can be neglected 
and one can use the equation ^48(1. We are interested in 
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the asymptotic behaviour of the polymer size moments 
M q (t) = (R q (t)) oc exp(A q t). In this case the value of 
a will be determined from the equation A a = 0. Inte- 
grating out the function £ p one can rewrite M q in the 
following form: 

M q = cxp [Dq 2 t - jqt] J dcj)Z q (q>, t) (54) 

Z q = (exp[ M (t)]5(0 - 0(t))) (55) 

where the angle brackets here stand for the averaging over 
the process 4>(t). The function Z q obeys the equation 



d t Z q 



2Ddl + sd 6 sin 2 , 



qs 



sin 2<f> 



(56) 



The only difference from the Fokker-Planck equation 1101 
only in the last term. We will follow the same procedure 
as in the section lTlI Bl Substituting Z q — U (0)5'(t, <f>) we 
obtain the imaginary time Schroedinger equation 

di* = -H q * (57) 

s 2 

H q = -1Dd\ + — sin 4 4> + (q - l)s sin cos 4> (58) 

This equation can not be solved in the case q ~ 1, how- 
ever for g > 1 the solution can be easily found. In this 
case the main exponential asymptotic at large times is 
Z q (t) —oc exp(— e(q)t), where e(q) is the ground state en- 
ergy. For g > 1 the main contribution to e(q) will be 
equal to the value of classical minima of the potential. 
After some simple algebra we obtain 



e(q) ~ -3 • 2- 5 / 3 (q 4 Ds< 



,1/3 



1 < g< e/D(59) 



e(q) ~ qs/2, q > s/D 



Finally we have 



A(q) = Dq 2 - 7 g - e(q) 



(60) 



(61) 



and the critical value a will depend on the dimensionless 
parameter 7 /s: 



81 7 3 

321^2' 

7 

a= D- 



7 -C s 
7 » s 



(62) 
(63) 



The last expression (|63|) coincides with the value of ex- 
ponent for pure isotropic chaotic velocity, because in the 
case 7 3> s large polymer size fluctuations are deter- 
mined by the rare fluctuations of the chaotic component, 
when the flow has a strong elongation component with 
the lyapunov exponent A > 7 for a long time. However 
the result 1|62JI shows that in the case 7 <C s the shear 
component can significantly broaden the tails of the poly- 
mer size PDF. This fact is nontrivial because the regular 
shear component itself can not lead to the exponential 
polymer elongation, and non-trivial exponent comes from 
the combined effect of chaotic and regular component. 



V. CONCLUSIONS 

In conclusion we would like to repeat the main re- 
sults of this paper. A polymer molecule injected in the 
chaotic flow with the strong mean shear component be- 
comes strongly stretched in the shear direction and expe- 
riences the a-periodic tumbling in the plane of the shear 
flow. It is shown, that stationary distribution of the poly- 
mer orientation angles has intermediate algebraic asymp- 
totics. The tails of the </>-angle distribution behave like 
</>~ 2 and correspond to the constant flux solution of the 
Fokker-Planck equation. Existence of these asymptotics 
has been reported in previous papers regarding the dy- 
namics of molecules and small particles in the shear flows 
[TsU^ . However the algebraic asymptotic of the 0-angle 
distribution has not been reported in any previous pa- 
pers to our knowledge. We would like to stress, that in 
contrary to the PDF of the 4> angle, the asymptotic be- 
haviour of the 0-angle distribution is not universal and 
depends on the statistical properties of the chaotic veloc- 
ity component. Therefore in our opinion measurement 
of the scaling exponent of the #-angle distribution is a 
perfect tool for the investigation of the statistical prop- 
erties of the chaotic velocity. Another quantity which 
can be measured by direct polymer observations is the 
PDF of the time periods between consequent polymer 
tumblings. The equilibrium direction of the polymer, 
which is oriented along the shear direction is semistable. 
This property leads to frequent polymer tumblings, so 
that the characteristic time between tumblings is of order 
r ~ r c = (Z?s 2 ) -1 / 3 . The situation would be different in 
the case of vesicles or non-spherical particles in the shear 
flow |27| , where the tumbling rate would be exponentially 
small. In the fourth part of this paper we make the pre- 
diction on the polymer size distribution. We show, that 
the existence of the strong shear component leads to the 
significant broadening of the size distribution compared 
to the isotropic case, which was considered in [||. This 
effect is to our opinion rather non-trivial, because the 
shear component itself can not lead to the exponential 
elongation of the polymer and the distribution broad- 
ening is the combined effect of the chaotic and regular 
velocity components. Finally, we would like to mention 
that all results of this paper were obtained under the 
assumption of isotropic and short-correlated chaotic ve- 
locity flow. While the first assumption, as it was dis- 
cussed throughout the paper is irrelevant for most of the 
results, the finite correlation time of the velocity flow can 
lead to the change of some quantitative results of this pa- 
per. However we believe (see [l4| for the more detailed 
discussion) that most qualitative predictions will remain 
valid if the velocity correlation is smaller or comparable 
to the characteristic time scale of the polymer dynamics 
t c . 
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